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ABSTRACT 

Penguin effects in the CP asymmetries of B® — > - , B Q d -> p ± 7rT and B Q d -> 
ct!j' = 7r~'~ are studied as function of the CKM unitarity triangle a. Despite a fairly 
small penguin amplitude, it leads to quite sizable uncertainties in the determination 
of sin(2a) from all but very large asymmetries. This effect is maximal for vanishing 
final state interaction phases, for which it can cause, for instance, an asymmetry 
of 40% if a = tt/2. 
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There are two good reasons for which B mesons provide a unique opportunity 
for testing the Cabibbo-Kobayashi-Maskawa (CKM) mechanism of CP violation 1 
. The CP asymmetries in certain decays, most notably decays to CP- eigenstates, 
are expected to be both large and theoretically clean. For instance, in B® — > ipKs 
the time-dependent asymmetry is predicted to oscillate with an amplitude given 
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directly by sin(2/3) , where (3 is one of the angles of the CKM unitarity triangle. 
This relation between a measured asymmetry and a pure CKM phase parameter 
follows from having essentially a single weak phase which contributes to the decay. 
In the case of B^ — > ipKs, where it is known that sin(2/3) > 0.08 1 , this relation is 
expected to hold within a 1% accuracy . 

Another decay mode which seems to be very promising is B® — > , in which 

the asymmetry is related to the angle a. In this case the theoretical situation is 
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somewhat less clean due to the contribution of "penguin" amplitudes , which 
may interfere with the dominant "tree" amplitude through their different weak 
phases. Here the ratio of penguin-to-tree amplitudes is roughly estimated to be at 
the level of (10 — 20)%. This estimate may lead one to conclude that an uncertainty 
at this level applies also to the relation between the measured asymmetry and 
sin(2o;). This by itself would not have spoiled the testing power of the asymmetry 
measurement. The purpose of this note is to critically elaborate on this question. 
We will show that, even with a relatively small penguin contribution and with 
small final-state interaction phases, penguin effects on the asymmetry may, in 
fact, be quite large for | sin(2ai)| < 0.7. Therefore, unless a very large asymmetry 
is measured, this would prohibit obtaining a useful value of a from the asymmetry 
measurement. 

The formalism of studying CP asymmetries in neutral B decays to CP eigen- 
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states in the presence of two interfering decay amplitudes was given in . For 
completeness, we write down the basic equations and study them for the case 
of B® — > 7r + 7r~. We denote the amplitude of B d — > it + -k~ by A and that of 
B d — > 7r+7r~ by A. Each of these amplitudes obtains contributions from "tree" 
and from "penguin" amplitudes: 

A = A T e iST e i<t,T + A P e i5p e^ p , 

A = A T e i5r e~ i<t)T + A P e iSp 'e~ i(t>p . (1) 

At,p are real, 4>t,p are CKM phases and 8t,p are strong interaction final state 
phases, all corresponding to the "tree" and "penguin" amplitudes, respectively. It 
should be mentioned at this point that 8t,p stand for soft final state interaction 
phases. We neglect a phase due to the absorptive part of the physical cc quark 
pair in the penguin diagram 6 . This phase is very small at the inclusive b — > uud 
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level , and is not expected to be considerably larger for exclusive modes such as 
7r + 7r~, where the absorptive part picks up contributions from a limited q 2 range 



The time-dependent CP asymmetry for a neutral B meson, created at t = as 

a pure B® and decaying at time t to n + 7v~ , when compared to the corresponding 
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decay rate of an initially B d , is : 

1 - |2| 2 2Im(2 e -2*/3) 

Asym{t) = ^- cos(Amt) -^-= sin(Amt) . (2) 



Am is the mass- difference of the two neutral B mesons. The phase 2f3 appears in 
the B® — B° d mixing amplitude in the standard CKM phase convention. The two 



terms in Eq.(2) describe two different kinds of CP violating phenomena. The first 
cos(Amt) term, which describes CP violation in the direct decay of B®, vanishes 
when only a single CKM phase contributes to the decay process. The coefficient of 



the well-known sin(Amt) term, which appears when the mixed B^ and B d decay to 



a comomn final state, is given by sin(2o;) when only one CKM phase contributes. 
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For Ap/Ap <C 1 one finds the following expressions for the two coefficients in 



Eq.(2): 



1-|-| 2 A P 

=— ~ -2-j— sin(0 T - <pp) sin(5 T - 5 P ) , 

1 + |4|2 A T 



2Im( A e -2iP) A 

y -^-= « sin(2a) + 2-f- sin(0 r - <f> P ) cos(2a) cos(5 r - 5 P ) . (3) 

1 + |A| 2 A T 

In order to evaluate the penguin effects one must know the three quantities (<pp — 
(ftp), Ap/Ap, (St — S P ). Only the first quantity can be studied theoretically in a 
reliable manner. In the standard CKM phase convention <f)p = phase(V* b V u< i) = 
7. The penguin amplitude, on the other hand, obtains contributions from three 
diagrams in which u, c, t quarks run in a loop. Denoting the three amplitudes, 
from which the CKM factors are omitted, by P u , P c , Pt, one notes that within 
0(ml/ml) one has P c ~ P u - It then follows from the unitarity of the CKM matrix 
that <pp ~ phase(V t * b V t d) = —(3, and therefore 

<pT — 0p~7 + /5 = 7r — a . (4) 



The ratio Ap/Ap cannot be calculated reliably at present. We will attempt to 
evaluate it in two different manners. To be conservative, we will try not to over- 
estimate it. The ratio of loop-induced processes b — > dqq (q = u,d, s,c),dg,dgg 
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to the tree process b — > uud was calculated perturbatively at the quark and gluon 
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level for the rates of inclusive charmless-strangeless decays . The penguin pro- 
cesses b — > dqq dominate the loop induced processes. The penguin-to-tree ra- 
tio of rates was found to decrease as a function of \y u bl^cb\ within the range 
0.07 < \V u b/V c b\ < 0.20, from a largest possible value of 0.4 to a lowest value of 
0.1, depending on m t and on CKM parameters. The quark process b — > duu is 
likely to be the dominant mechanism for B® — > . Its rate is smaller by about 
0.3 relative to all b — > dqq processes, which leads to a penguin-to-tree ratio of 
amplitudes decreasing from 0.35 to 0.15 in the above range of |V^&/Vrf,|. This may 
be a slight overestimate, since in this calculation the penguin rate was maximized 
over the entire acceptable CKM parameter space. 

A somewhat different approach which leads to a similar estimate of Ap/At is 
based on calculating the low energy effective Hamiltonian for b — > duu in the lead- 
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in log approximation . One finds that the penguin operators du^bLjULjl 'pUu, 
dLil^bijURj'j^URi, are color indices and L,R are left and right projections) 
appear with coefficients 0.026, 0.033, respectively, and are multiplied by the CKM 
factor V t bV t * d . On the other hand, the relevant tree operator has a slightly en- 
hanced coefficient (= 1.11) and is multiplied by V u bV* d . Adding up the penguin 
coefficients and allowing the ratio I^V^/l^V^I to vary between the values of 1 
and 5 1 , one finds at the quark level Ap/A T ~ 0.05 — 0.27. Note that, since the 
operator coefficients are all positive, Ap/Ap (which does not include final state 
phases) is positive too. 

A very rough and oversimplified approximation, which represents the above 
two results, somewhat on the low side, can be obtained by simply using the CKM 
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factors and the QCD factor related to the single t-quark penguin diagram: 



A P \V*V td \ 



as In A) ~ 0.04 - 0.20 . (5) 



A T \V*V ud \l2n W 



The ratio of the CKM factors is the ratio of the lengths of two sides of the CKM 
unitarity triangle which form the angle a. As noted above, this ratio lies between 
the values of 1 and 5 1 . Recent preliminary data 10 , which seem to indicate that 
\Vub/Vcb\ is only about 0.06 or even smaller, favor a large ratio. For the QCD factor 
(a s /127r) ln(m^/m^) we took the value 0.04, using a s (ml) « 0.2. This value would 
be larger if a s were to be taken at (m&/2) 2 . 

A large uncertainty is involved in calculating the tree and the penguin operator 
matrix elements between the and the ir + ir~ states 11 . Certain hadronic models 
seem to indicate that penguin operator matrix elements may be enhanced due to 
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their special chiral structure . Since none of the existing methods of calculating 
hadronic matrix elements is very reliable for our case, we will make the most 
simplified assumption that the ratio of these matrix elements is one, and will 
thus use Eq.(5) as a crude approximation. This assumption has not yet been 
tested experimentally even in an indirect way, that is, by comparing tree-dominated 
to penguin-dominated processes. We feel that, since the simplified relation (5) 
somewhat underestimates the ratio calculated at the quark and gluon level, it 
allows a certain amount of penguin matrix element suppression, and does not 
overestimate the ratio of matrix elements. We note again that this ratio is likely 
to be on the high side of (5) if IV^/V^I is near its present lower limit value of 0.06. 

The soft final state interaction phase difference, 5t — Sp, is basically uncalcu- 



lable. Denoting this phase- difference by 8, one finds from (2)-(4): 

A P 

Asym(t) ?s — 2— — sin5sino;cos(Amt) 

Ap 

(6) 

Ap 

— [sin(2a) + 2— — cos<5cos(2o;) sina]sin(Amt) . 
Ap 

We note that the cos(Amt) term and the sin(Amt) term have a different and 
complementary 5-dependence. The first term, which describes CP violation in 
the direct decay, behaves like sin 5, while the correction to the mixing-induced 
asymmetry is proportional to cos 5. Thus, as function of 5, the smaller the direct 
CP violation cos(Amt) term, the larger becomes the penguin correction to sin(2o;), 
and vice versa. In particular, when 5 = 0, the cos(Amt) term vanishes, whereas 
the correction to the sin(2o;) coefficient becomes maximal. 

A heuristic argument for factorization of tree amplitudes in certain two body 
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B decays implies that 5p is negligible and that perhaps also Sp is small. // 

5 were small, then the cos(Amt) term may be too small to be observed and the 
time- dependent asymmetry measurement would not provide evidence for a penguin 
contribution. Still, in this case the penguin amplitude effect on the coefficient of 
the sin(Amt) asymmetry becomes maximal and may be large. This is the danger 
of penguin amplitudes. 

The crucial point is that, whereas one would naively expect that the penguin 
amplitude modifies sin(2o;) in a multiplicative manner, the correction is in fact 
an additive one and involves a factor of two from the interference with the tree 
amplitude. This means that with e.g. Ap/Ap = 0.2, the correction to sin(2a) can 
be as large as ±0.4 and is not merely a relative 20% correction. To be quantitative, 
let us assume that 5 is negligibly small, and study the consequences of (6) on the 
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determination of sin(2a) from an asymmetry measurement. Fig. 1 shows the 
coefficient of the — sin( Amt) term as function of the actual value of sin(2o;) for 
45° < a < 135°. The range bounded by the two solid lines describes this coefficient 
for 5 = and for Ap/Ap in the range (5). The straight dashed line gives the 
corresponding relation in the absence of a penguin contribution. The maximal 
deviation from the straight line is given by 2Ap/Ap- We note that points with 
largest deviations from this line correspond to the largest value of Ap/Ap and thus 
to the smallest values of |V^&/T<^&|. The danger of the penguin amplitude is best 
demonstrated for a = ir/2, where its effect on the asymmetry is maximal. For this 
case, an asymmetry as large as 0.4 can possibly be measured, although sin(2o;) = 0. 
Such a substantial CP asymmetry measurement would be an important observation 
by itself, however it could not be related to the angle a. The deviation of the 
asymmetry from sin(2a) decreases gradually, as one moves away from a = n/2. It 
is still at a level of 30% at a = 65°, 115° (sin(2a) = +0.77, -0.77), where one 
expects large asymmetries. The effects become much smaller outside the range 
45° < a < 135° plotted in Fig. 1. Unfortunately, the asymmetry measurement, 
which is related to the value of sin(2o;), cannot distinguish between angles which lie 
outside and inside this range. Note that as far as the CP asymmetry measurement 
is concerned, the corrections are potentially large for all but very large asymmetries. 
The range 65° < a < 115°, where the corrections are larger than 30% 14 , is 
presently allowed 1 . It is possible that future theoretical and experimental progress 
in quantities such as fp (the B decay constant), (V^/V^l and m t will rule out this 
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range 

If one takes the very conservative viewpoint that the sign of the penguin ampli- 
tude relative to the tree amplitude is unknown, then the uncertainty in determining 
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a from the asymmetry becomes twice as large. Reversing the sign of Ap/At cor- 
responds to flipping the allowed range for the asymmetry coefficient to the other 
side of the straight dashed line. It would therefore be useful to at least theoreti- 
cally determine the sign corresponding to 5 < ir/2. Since the penguin correction 
to sin(2a;) is proportional to cos 5, it becomes substantially smaller than in Fig. 1 
only for large values of the final state phase difference. In this case one expects to 
observe also the cos(Amt) term in the time-dependent asymmetry. This depends, 
of course, on the value of a and on the sensitivity of the experiment. With no 
observation of such a term, and without a theory of final state interaction phases, 
one would have to assume the worst of all cases, namely 5 ~ (or even 8 ~ ir, if 
the sign of the penguin amplitude is undetermined). 

Penguin contributions appear also in B d — > p ± 7T =F and in B d — > afir^. The 
effects on a determination of a from the asymmetry of the related time-dependent 
rates is expected to be as large as in B d — > it + it~. We wish to briefly demonstate 
this effect in the decays to pir. The general formalism dealing with this kind of 
final states (which, although being non-CP-eigenstates, are common to B d and B d 
decays), was described in 16 . The essential difference with respect to B d — > 7r + 7r~ 
is that here the tree and penguin amplitudes for B d — > p + n~ are not the same as 
those for B d — > p~n + . In both cases the corresponding amplitudes carry the same 
CKM phases as in B^ — > it + it~, 4>t = 1 and <pp = —f3, respectively: 

A f = A(B° d -> P+7T-) = A T e iST e i(f>T + A P e iSp e^ p , 

Aj = A{B° d -> p-7i+) = A T e f5T e ict,T + A P e fSp e i<t>p . (7) 
The corresponding amplitudes for the charge-conjugated processes, Aj = A(B d — > 

9 



p 7r + ), Af = A(B° d — > p + 7r ), are obtained simply by changing the sign of the 
weak phases. 

One can measure four different time-dependent decay rates, for cases in which 
initially B d (B^ d ) decay to p ± 7r =F . If Ap could be neglected then these four rates 
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would be, in principle, sufficient for a determination of a .As in decays to CP 
eigenstates, this method involves a measurement of the coefficient of the sin(Amt) 
term in the time-dependent rate. For an initially B d decaying to p + n~ this coeffi- 
cient is given by 

Im(4^e _2i/J ) = 4 1 sin(2a + 5 T - 6 T ) ■ (8) 
A f Ap 

If all final state phases were negligible, this coefficient would determine sin(2ai). 
(Otherwise, the phase difference 5p — 8p can be determined separately from the 
four rates). In the presence of the penguin amplitudes, this coefficient becomes, to 
lowest order in Ap/Ap and Ap/Ap, and for negligible final state phases, 

l m (4^ e - 2i/3 ) = 4 I [sm(2a) + 4^sin(3a) - =^ sin a] . (9) 
y Af Ap Ap Ap 

The correction to sin(2o;) can be estimated in a way similar to the correction 
in B® — > 7r + 7r~. Within our approximation, Eq.(5) applies to both Ap/Ap and 
Ap/Ap. In fact, if one takes Ap/Ap = Ap/Ap (which should hold only for CP- 
eigentates) the correction term obtains exactly the form of the correction term of 
Eq.(6). In general, when these ratios are in the range (5), the effect of the penguin 
amplitude on determining sin(2a) are expected to be as large as in B® — > . 

In summary, we have shown that relatively small penguin amplitudes may pro- 
hibit a useful determination of sin(2a) from the CP asymmetries of B® — > -k + it~ , 
— > af-K^ 1 . Asymmetries as large as 0.4 may be measured even 
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when sin(2o;) = 0. The uncertainty becomes small only for very large asymmeries. 
It would decrease if (V^/V^l were found to be on the high side of the presently 
allowed range. The penguin complication may be avoided to a large degree if future 
studies of the CKM matrix exclude the range 65° < a < 115° in which the correc- 
tions are large. Our analysis was based primarily on the estimate (5) and on the 
observation that the correction to sin(2o;) in (6) is additive rather than multiplica- 
tive, and becomes maximal when 5^0. We assumed that the hadronic matrix 
element of the penguin operator is neither dynamically suppressed nor enhanced 
relative to the tree amplitude. It goes without saying that this issue deserves se- 
rious studies, both theoretical and experimental. The mere determination of the 
sign of the penguin amplitude would be useful. 

One way to overcome this potential difficulty is to measure in addition to the 

asymmetry in B® — > also the rates of — > 7r°7r°, B + — > ir + n°. This isospin- 
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based method can provide a way to eliminate the penguin contribution altogether 
and to experimentally determine its magnitude, provided that the integrated rate 
into two neutral pions is measurable. A similar isospin analysis for the p7r modes 
is unlikely to work in practice due to the too many amplitudes involved and to 
certain ambiguities which appear in the analysis 19 . 

I wish to thank J. Bjorken and H. Quinn for useful discussions. 
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FIGURE CAPTION 

FIG.l. Coefficient of — sin(Amt) in the asymmetry of — > it + it~ as function of 
sin(2o;) for 45° < a < 135°. We take 5 = 0. Area between solid lines 
corresponds to Ap/At = 0.04 — 0.20; dashed line corresponds to the absence 
of a penguin amplitude. 
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